Agglomeration of microparticles in complex plasmas 
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Agglomeration of highly charged microparticles was observed and studied in complex plasma 
experiments carried out in a capacitively coupled rf discharge. The agglomeration was caused by 
strong dust density waves triggered in a particle cloud by decreasing neutral gas pressure. Using a 
high-speed camera during this unstable regime, it was possible to resolve the motion of individual 
microparticles and to show that the relative velocities of some particles were sufficiently high to 
overcome the mutual Coulomb repulsion and hence to result in agglomeration. After stabilising the 
cloud again through the increase of the pressure, we were able to observe the aggregates directly with 
a long-distance microscope. We show that the agglomeration rate deduced from our experiments is 
in good agreement with theoretical estimates. In addition, we briefly discuss the mechanisms that 
can provide binding of highly charged microparticles in a plasma. 
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I. INTRODUCTION 



Dust agglomeration, also known as dust aggregation 
or dust coagulation, is a process of sticking dust particles 
together after their collision [l|, [2| . This process occurs 
in very different environments and plays an important 
role in colloid chemistry Q , astrophysics [J] , atmospheric 
physics [5|, fusion physics Q as well as plasma physics 
[7|. Nano/microparticles trapped in a low-temperature 
plasma are normally charged negatively due to high elec- 
tron mobility [8|, which results in the mutual repulsion 
between them. Agglomeration of nanoparticles in plas- 
mas - the major mechanism operating during the so- 
called "rapid growth" (RG) phase Q - has been observed 
in numerous experiments jlOllllj . At the RG stage, most 
of nanoparticles are neutral, some carry several electron 
charges, even fewer are charged positively (Tl|, Qjl], and 
therefore the charge discreteness plays very important 
role. The agglomeration criteria for such particles were 
derived theoretically, based on either the polarization- 
induced ion flux asymmetry model [13| or the shadow 
attraction due to neutral/plasma particle bombardment 
[L4J . The agglomeration of nanoparticles in a plasma 
has been reported for sizes up to several hundreds of nm 
(la Ha ] . As the size of particles immersed in a plasma 
increases, their charges become more and more negative 
(the charge is linearly proportional to the particle size, 
for a one micron particle it ranges from several hundreds 
to several thousands of e, depending on discharge condi- 
tions [3, H3, [ill)- Therefore, agglomeration of micropar- 
ticles in a plasma has been believed so far to be practi- 
cally impossible due to strong Coulomb repulsion. This 
process was primarily studied (both experimentally and 
theoretically) in a neutral environment 19, 20], showing 
that positive and negative charges carried by individual 
particles (in overall neutral ensemble) can enormously 
accelerate the agglomeration. 

In this article we report on a series of experiments 
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where dust-acoustic (density) waves trigger agglomera- 
tion of microparticles. We show that the agglomeration 
can occurs under typical discharge conditions, for parti- 
cles of different sizes and materials, and the aggregates 
formed in a plasma can be observed in situ, using a long- 
distance microscope. 



II. EXPERIMENTS AND OBSERVATIONS 

The experiments were performed on the ground in the 
PK-3 Plus chamber [2l| . which was originally designed 
for microgravity experiments onboard the International 
Space Station (ISS). The chamber (which consists of two 
parallel electrodes driven in push-pull mode at a fre- 
quency of 13.56 MHz) was placed horizontally and sin- 
gle microparticles (no agglomerates) were injected into 
the discharge from a special dispenser. A vertical laser 
sheet (perpendicular to the electrodes) with a thickness 
of about 100 /im was used to illuminate the particles in 
the plasma, the scattered light was recorded at 90° by a 
high-speed CCD camera. A detailed description of the 
PK-3 Plus setup can be found in |2l|. 

Microparticles were first injected into argon plasma 
(typically maintained at a rf peak-to-peak voltage of 
~ 25 V and gas pressure of ~ 50 Pa), to form a sta- 
ble cloud above the bottom electrode. Then, in order to 
accelerate particles, dust-acoustic (density) waves were 
triggered in the cloud by decreasing gas pressure well 
below the instability threshold [18| (down to 10-20 Pa). 
This unstable regime was maintained for several minutes, 
and then the gas pressure was restored to the initial (sta- 
ble) level. Such experiments were performed with parti- 
cles of different sizes and materials. 

We compared profiles of the dust clouds before and af- 
ter the wave excitation, as illustrated in Fig. [TJ After the 
excitation, we observed aggregates levitated at the bot- 
tom of the cloud (they can be identified as brighter spots 
having larger distances to their neighbors). For Si02 
particles of 1.5 /im diameter (first row), multiple layers 
of aggregates (apparently of different sizes) were seen af- 
ter the unstable regime. The aggregates were separated 
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FIG. 1: Particle cloud before the dust density waves were triggered (first column) and after the waves were stopped (second 
column). Shown are results of the experiments with different particles: The first row represents Si02 particles of 1.5 fim 
diameter, the second row - MF particles of 3.4 /mi diameter, the third row - MF particles of 6.8 /jm diameter, and the fourth 
row - PMMA particles of 3 /im diameter. 



from single particles by a narrow but well-seen gap, while 
before the wave excitation there were only a few big- 
ger particles near the lower edge. The similar situation 
was also observed for melamine-formaldehyde (MF) par- 
ticles of 3.4 /jm and 6.8 /im diameter (second and third 
row, respectively) as well as for polymethyl-methacrylate 
(PMMA) particles of 3 /im diameter (fourth row). 

In order to take a close look at the structure of the ag- 
gregates we used a long-distance microscope (with a dif- 
fusive light source from the opposite side of the chamber), 
which allowed us to perform in situ observations. Due to 
limited resolution of the microscope, we were only able 
to analyze the experiments with sufficiently large parti- 
cles. The results are presented in Fig. [2] In the experi- 
ments with MF, PMMA, and "magnetic polystyrene" 22 1 
(MPS) particles of about 3 /im diameter (first, second, 
and third row, respectively) the aggregates were usually 



composed of several particles and had rather irregular 
(often elongated) structure. For larger MF particles of 
6.8 /im diameter (fourth row) the aggregates were typi- 
cally formed of two-three particles. It is noteworthy that 
for silver-coated polystyrene (PS) particles of 5 /im diam- 
eter (fifth row) the agglomeration was also possible and 
in fact very effective, resulting in the formation of rela- 
tively large elongated aggregates. This strongly suggests 
that the agglomeration of microparticles in a plasma is 
a quite generic phenomenon occurring both for dielectric 
and conductive particles. 



III. ANALYSIS 

The agglomeration in a plasma can only be achieved 
if heavily charged particles acquire the kinetic energy 




FIG. 2: Single micro-particles before the waves were triggered 
(first column) and examples of aggregates after the waves were 
stopped (second and third column), all observed with a long- 
distance microscope. Experiments were performed with dif- 
ferent particles, the first row shows MF particles of 3.4 /xm 
diameter, the second row - PMMA particles of 3 /im diam- 
eter, the third row - MPS particles of 3 /jm diameter, the 
fourth row - MF particles of 6.8 /im diameter, and the fifth 
row - silver-coated PS particles of 5 fj,m diameter. 



which is sufficiently high to overcome the interparticlc 
repulsion. In our experiment, the particles were acceler- 
ated by the dust-acoustic waves (which were presumably 
triggered by the ion streaming instability [8j). In or- 
der to track the individual particles inside the wave, we 
performed a series of experiments with MF particles of 
2.5 /xm diameter and recorded their motion during two 
seconds. We used a high-speed camera at 2000 frames 
per second (fps) with full exposure, recording the field of 
view of 3 mmx3 mm (at a spatial resolution of 210 pixel 
per mm). The waves propagated mostly vertically, with 
an average wave period of about 60 ms (at a gas pressure 



of 20 Pa). The waves were highly nonlinear, with a mi- 
nor fraction of particles moving in one direction at rather 
high velocities, while most of them were moving slowly 
in the opposite direction. To visualize the particle tra- 
jectories, five consecutive video frames were color-coded 
and overlapped, as shown in Fig. |3Ji. 

In order to analyze the agglomeration process in de- 
tail, we first have to evaluate the parameters that deter- 
mine the interparticle interaction. Obviously, the most 
critical parameter for such analysis is the charge of mi- 
croparticles - if it is too high the coagulation is inhib- 
ited. A simple estimate of the charge from the orbit- 
motion-limited (OML) theory is just not possible for our 
conditions, since the plasma parameters in dense parti- 
cle clouds are poorly known. As an alternative, we were 
able to derive the charge from rarely occurring scattering 
collisions of two particles. We supposed that the ob- 
served collision events arc representative, in the sense 
that charges of the colliding particles are close to the 
mean particle charge in the whole cloud. An example of 
the collision event is shown in Fig. [3]>f. For the anal- 
ysis we assumed the Debye-Hiickel interaction potential 
which depends on two unknown parameters: the charge 




FIG. 3: Trajectories of individual microparticles in the wave, 
(a) The trajectories are depicted by superimposing five con- 
secutive video frames, color-coded from dark blue to red. (b)- 
(f) Example of a two-particle scattering collision [magnified 
rectangular box from (a)], the scattering is shown as a se- 
quence of five consecutive frames. The results are for MF 
particles of 2.5 fim diameter. 



Q and the screening length A. These parameters were 
derived by simulating the pair collisions and fitting the 
resulting trajectories to the observed ones. We took into 
account neutral gas friction |23J and used the particle 
velocities and coordinates measured before the collision 
as the initial conditions. Since the velocities and coordi- 
nates were naturally measured with some errors, the fit- 
tings of Q and A was performed by varying the initial con- 
ditions within the estimated error bars. The best fit was 
then defined by maximizing the correlation between the 
simulated and measured trajectories, which is illustrated 
in Fig. |U We found that the fit is practically insensible 
to the value of A, apparently because the relative veloci- 
ties of colliding particles were so high that the scattering 
mostly occurred at distances smaller than the screening 
length. This observation significantly simplified the sub- 
sequent analysis and allowed us to directly employ the 
Rutherford's scattering theory [IS, l2J|- Given the fact 
that the typical collision event is much shorter than the 
gas damping timescale, the charge can be obtained from 

the formula Q = — i/2.E&tan|0, where E = jmv 2 is 

the (initial) kinetic energy of the relative particle motion 
(where m is the mass of a single particle and v is the 
relative velocity), 9 is the scattering angle, and b is the 
impact parameter for colliding particles. This yields the 
charge Q = —(1650 ± 170) e, where the uncertainty is 
given by the standard deviation averaged over ten colli- 
sion events. However, one should keep in mind that the 
actual uncertainty is somewhat larger, because the colli- 
sion plane (in the center-of-mass reference frame) is never 
perfectly parallel to the illuminating laser sheet (which 
has a finite thickness), as one can see from Fig. |4] 

Two colliding particles can touch each other only if 
their relative velocity exceeds a critical (threshold) value 
determined by the Coulomb repulsion. Using the same 
assumptions as above, we readily obtain the critical ve- 
locity, 
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which has the following scaling on the particle size: v cv oc 
a" 1 (here we also assume Q cc a). For the experiment 
shown in Fig. [3] we get f cr ~ 273 mm/s. In order to 
evaluate the probability for particles in the cloud to reach 
the critical velocity, we tracked the individual particles 
manually and estimated their velocities from the lengths 
of the trajectories at each video frame. From this analysis 
we deduced the relative velocity distribution (normalized 
to unity), /rci(w), which is shown in Fig. [5] for different 
phases 4> of the wave. This plot clearly demonstrates that 
for certain <f> there is a small fraction of particles with the 
relative velocities exceeding v cr . 

Using these results, we estimated the agglomeration 
rate for pair collisions (per unit volume, assuming that 




FIG. 4: Comparison of a two-particle scattering collision with 
the simulations. The simulated trajectories (marked by dia- 
monds and squares) are overlapped with the observed particle 
trajectories (from the experiment shown in Fig. [3]). One of 
the trajectories (diamonds) becomes noticeably dimmer at 
one end, suggesting that the scattering plane (in the center- 
of-mass reference frame) is not parallel to the illuminating 
laser sheet. 



the sticking probability upon the collision is unity), 

1 f T f°° 

R= — dtn 2 dv o-y/rei, 

1 Jo Ja 
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where n(t) is the (time-dependent) particle number den- 
sity in the cloud, T is the wave period, and the collision 
cross section, 
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is determined by the critical velocity given by Eq. ([T}. 
The calculations yield R ~ 0.8 mm _3 s _1 for the ex- 
periment illustrated in Figs [3][5] In order to compare 
these estimates with the experiment, we selected a part of 
the visible cloud (assuming the horizontal homogeneity) 
which represent a slice of the thickness ~ 0.1 mm (about 
the thickness of the illuminating laser sheet), width of 
3 mm, and height of 5.8 mm (full height of the cloud). 
The number of (double) aggregates predicted by Eq. © 
to appear in this volume after five minutes of the unsta- 
ble regime is ~ 450, which is 5-6 times larger than the 
number of aggregates (82) we actually observed. Taking 
into account that the sticking probability of the collid- 
ing particles can be substantially smaller than unity and 
also that the aggregates in the considered experiment 
consisted (on average) of more than two particles, we 
conclude that the theory provides reasonable agreement 
with the experiment. 



IV. DISCUSSION AND CONCLUSION 

Let us discuss the mechanisms that can effectively bind 
particles in a plasma. If we assume a homogeneous charge 
distribution over the spherical particle surfaces, the re- 
pulsive Coulomb force Fq between them is equivalent to 
that between two point charges Q located at the centers 
of the spheres. Therefore, the force at the contact, 
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does not depend on a. On the other hand, the charge at 
one particle induces a dipole moment at another particle. 
This results in the attractive force which depends on the 
permittivity of the particle material and also scales as 
oc (Q/a) 2 , thus tending to balance Fq. 

However, there are (at least) two other very efficient 
attraction mechanisms that can overcome the repulsive 
force (J3J) - One of them is due to the van der Waals inter- 
actions between individual molecules of microparticles. 
For two identical spherical particles, Hamaker |25| gives 
the following expression of the attractive force: 



-FvdW = 



4« 
12<S 2 ' 



(4) 



where Ah is the material-dependent Hamaker constant 
and S(<^ a) is the effective distance between surfaces 
of two particle at the contact 0, (25|. (If the polariza- 
tion of spheres is neglected then Ah = ir 2 p 2 a, where 
p is the number density of molecules and a is the cor- 
responding London- van der Waals constant.) By com- 
paring Eqs ((3|) and (j4]) we see that the van der Waals 
attraction is stronger than the Coulomb repulsion when 
Ah > Q 2 5 2 /a 3 ((x a -1 ). Let us consider an aggregate 
composed of two Si02 particles (An = 6.5 x 1CT 20 J) of 
the diameter 1.5 pm (the smallest particles used in the 
experiments, see the first row of Fig. [1}. Assuming that 
the particles carry the charge Q = —1000 e each and are 
separation by S ~ 3 A at the contact, we obtain that 
-Fvdw is two orders of magnitude stronger than Fc . Note 
that for metal particles (such as Au, Ag, Cu) Ah is typi- 
cally in the range (25 — 40) x 10 -20 J |26J |. and therefore 
the binding force for them should be even stronger. One 
should also point out that the particle coating (see, e.g., 
the fifth row of Fig. [2]) may affect van der Waals attrac- 
tion, provided the coating thickness is comparable to or 
exceeds 5 and the a constant of the coating material is 
different from that of the particle material. 

The second mechanism of attraction can be associated 
with the discreteness of the charge distribution over the 
particle surface (here, we do not discuss the physics of 
the electron states at dielectric surfaces) . The "mean dis- 
tance" between the elec trons at the surface, which can be 
estimated as A e ~ ay / 4ire/\Q\ oc a 1 / 2 , is typically of a 
few dozens of nm. Therefore, when particles get in con- 
tact their interaction should strongly depend on whether 
there are individual electrons in a close proximity of the 
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FIG. 5: Relative velocity distribution of microparticles, 
/rei(f), for different phases <f> of the wave period. The dis- 
tribution is normalized to unity, the dash line corresponds to 
the critical velocity v C i(= 273 mm/s). 



contact spot. For instance, if there is one electron at 
each particle surface near the contact spot, this will ex- 
ert an additional repulsive force which can be as large 
as ~ e 2 /5 2 . However, if an electron is only attached to 
one of the particles, this results in the additional attrac- 
tion between them due to the local polarization of the 
particle material. Simple calculations (representing each 
particle as an infinite half-space filled with a material of 
high complex permittivity) yield the following attractive 
force due to the "charge discreteness" : 
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By comparing Eq. ([5]) with Eq. © we obtain that the 
charge discreteness force should overcome the Coulomb 
repulsion when [a/5) 2 > (Q/e) 2 , and this condition does 
not depend on a. For Si0 2 particles considered above 
(to illustrate the van der Waals attraction), Faisc is an 
order of magnitude stronger than Fq and hence can also 
contribute to the binding. Note that -Fdi sc is really of 
short-range and can only be important when particles 
get very close to each other (S < A e ) - otherwise the 
interaction is described by usual "mean-field" formulas. 
Also, iodise must be strongly dependent on the particle 
(surface) conductivity and presumably be much larger 
for dielectric particles than for conductive (or coated) 



One should make a common remark concerning both 
binding mechanisms discussed above. In our calculations 
we used rather crude assumption that the charges of in- 
dividual particles in a plasma and charges of particles 
bound in a cluster are the same. This assumption is 
quite justified as the "first approach" , since the estimated 
binding forces turn out to be much stronger than the 
Coulomb repulsion. However, for a more careful analysis 
of the problem one should take into account the mutual 
influence of the neighboring particles on their charges 
(see, e.g., [27], [28|]). In fact, the binding force should 
strongly depend not only on the magnitude of the parti- 
cle charge, but also on the charge distribution over the 
particle surface. This very complex issue undoubtedly 
requires further careful analysis. 

To summarize, by triggering dust-acoustic waves in a 
cloud of microparticles suspended in the pre-sheath re- 
gion of a low pressure rf discharge (which is referred to 
as the unstable regime), we were able to accelerate the 
particles to velocities which were sufficiently high to over- 
come their mutual Coulomb repulsion. This resulted in 
the agglomeration - the formation of clusters consisting 
of two, three, or more particles - yet this effect was ob- 
served for particles of different sizes and materials. By 
tracking individual particles and resolving their scatter- 
ing collisions with a high-speed camera, we calculated the 
mean particle charge and reconstructed the relative ve- 
locity distribution function. This allowed us to estimate 



the coagulation rate during the unstable regime and com- 
pare it with the experimental observations. Additionally, 
by using a long-distance microscope we were able to ob- 
serve the aggregates in situ, levitated in a plasma after 
the wave excitation. 

The results presented in this paper unambiguously 
demonstrate that microparticles, even exposed in a 
plasma environment, can form stable clusters due to 
strong short-range binding forces. We believe that these 
results can be of great importance to various fields rang- 
ing from plasma processing (where this effect can be uti- 
lized, e.g., to remove dust contamination) to astrophysics 
(where the effect can have a strong impact on the coag- 
ulation of cosmic dust). 
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